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Abstract 

In the present work we investigate the existence and stability prop- 
erties of q-balls which consist of a couple of scalar fields, forming an 
SU{2) doublet in a Lagrangian with a global SU(2) symmetry. We 
find that these spinors can form a localized and stable field configura- 
tion, if they rotate in their internal SU(2) space. We find the energy 
and charge of the soliton in both thin and thick-wall approximation 
and we prove its stability against decaying to free particles. We also 
find the asymptotic forms of the scalar and gauge field and the energy 
and charge of the configuration when the SU(2) symmetry is local. 
The only assumption is the smallness of the coupling constant g. Us- 
ing numerical methods we prove the stability of the q-ball in the local 
case. 
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1 Introduction 



Q-balls are non-topological solitons appearing in Lagrangians of scalar fields 
with a global U(l) or SO (2) symmetry. Though they appeared as a rather 
mathematical object £Q, through a series of papers they revealed their phys- 
ical properties [3 El El- Q-balls can be observed at, or near, the minimum of 



the yU/(\$\) quantity, where U is the potential and $ is the scalar field. At 
the above minimum one can prove that the energy is minimized with respect 
to the charge. The minimum of this quantity should be less than the particle 
mass, in order the soliton to be stable. Much more interest was concentrated 
on the subject when their possible existence in supersymmetric extensions 
of the Standard Model was considered 00 HI IE], where the £7(1) charge is 
for example the baryon number of the supersymmetric partners of baryons. 
Flat U(l) directions in the potential of such extensions can offer a plausible 
explanation to the problem of baryogenesis ITUj). 

Another type of q-ball type solutions has been studied, namely the so- 
called non-abelian q-balls [TT) IT2*] where we have scalar fields in a Lagrangian 
with a global 577(3) or 5*0(3) symmetry usually and the fields belong to the 
adjoint representation of the symmetry group. Such particles may correspond 
to the sgauginos in N = 2 supersymmetric theories. Also, q-balls with a local 
U(l) symmetry have been studied in [T3*] . 

In a recent series of papers [TJJ EH], a great number of SU{2) flat 
directions in the superpotential of the Minimal Supersymmetric Standard 
Model and a rich phenomenology concerning them have been investigated. 
So the plausible question is to ask if a new type of q-balls can appear in 
such directions, enriching the particle phenomenology. In the present work 
we try to do the following three things: a) We present the general formalism 
concerning q-balls when the fields transform according to the doublet repre- 
sentation of the symmetry group because this case has not been studied and 
may turn out to be useful in the search of q-solitons composed of supersym- 
metric partners of the usual leptons. b)We present the formalism of gauged 
q-balls which describe scalar particles with "weak" interactions. c)We fully 
investigate the properties of the above solitons both analytically and numer- 
ically in the so-called thin- and thick-wall approximation. Our main interest 
is to verify their stability when varying the gauge coupling, the frequency 
with which the soliton field rotates within its internal SU (2) space, or other 
phenomenological parameters. 
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2 Global SU(2) symmetry. 
Thin-wall approximation 

We regard a Lagrangian 

L = drfdPtp-U{$0) , (1) 

characterized by a global SU{2) symmetry concerning the couple of scalar 
fields: 

The energy is provided by the equation 

E 



J d 3 x(V^V0 + $<p + U) . (3) 

The currents are three as the generators of the algebra considered. These 
currents are 

£ = %' _\ T , ) ( |) , (4) 

where r a are the Pauli matrices: 

Tl = ( °i I ) T2 = ( °t ~o ) Ts = ( o -°i ) • (5) 

We choose 

(^(f,t) = exp (^r a e a tj B(r) . (6) 

In the so-called thin-wall approximation we choose a(r) a step-function, with 
a constant value within a volume V and zero outside. We also choose a 
spherically symmetric configuration, i.e. a(r) = a(r) so as to minimize the 
contribution of the surface effects to the total energy. So the energy can be 
written 

E = $0V + UV (7) 
and using the chosen field configuration: 

E= ^(el + el + eD^aV + UV . (8) 
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We can now define the charges in the above thin-wall approximation 
where we take 

Qi = \e^BV . (9) 



where the definitions: 



and 



Qa = I d 3 xj° a (10) 



Q = (Q 1 ,Q 2 ,Q 3 ) (11) 



hold. We now want to write energy as a function of the charges. It is a 
matter of simple algebra to show that: 

Regarding energy as a function of the volume and minimizing it with respect 
to the volume, we take that the minimum is at: 



v =«mu (13) 



and the energy takes the form: 



The energy of the free particles carrying charge y/ Q 2 is 2m • yj Q 2 . (It 
is a matter of simple algebra to prove this. The physical interpretation is 
obvious because we now have two particles, each with mass m.) So if the 
quantity (U/a^a) 1 ^ 2 is less than the mass of the free particles then the energy 
of the ()14|) field configuration is less than the energy of the free waves and 
the soliton is preserved stable against decaying into free particles. We know 
that a/Q 2 ", quantity is both conserved, as every single current is conserved, 
and SU(2) invariant. We now investigate a potential of fourth power in 
\cp\ = (<^<y2) 1//2 , namely: 

TT 2| |2 2a, ,3 p. ,4 /-, K \ 

U = m \(p\ - —\tp\ + -\<p\ . (15) 

This is the more general potential of fourth order in the fields. The theory 
is non-renormalizable but this is not a serious problem as we can regard the 
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Lagrangian as a phenomenological one. This potential can be proved very 
useful in a supersymmetric theory. Making the rescallings \(p\ — > (m 2 /a)\ip\, 
x — > x/m and defining B = f3m 2 /a we take 

TT I |2 2. .3 5 . .4 , . 

^ = M - + -77M ■ (16) 

Minimizing the (4U/\(p\ ) quantity with respect to the field and calling it 
\fe^ we can take that the minimum is at 



Then: 



\<p\=2/3B. (17) 



el = 2^1-2/95 , (18) 



Q = (1/2) (2/35) Ve , (19) 



E = (2/35) V ^el y/l -2/95 , (20) 
where e = (£i,£ 2 ,£3)- It is easy then to verify that E = a/s^a/Q^. 

3 Local SU(2) case. General setting 

The Lagrangian for the case with a local 577(2) symmetry is written 

L = ] (V + ig^^j fa - w-^^j 0-U(0*0) - \f^f, v , (21) 

where: 

F"" = d"A" - d v A" + gA 11 x A v . (22) 

and g is the coupling constant. 

We now regard: a) Spherically symmetric fields in order to minimize the 
contribution to the total energy of the so-called surface effects, b) Static 
gauge fields. c)Only the zero component of the gauge field differing from 
zero, the other components being equal to zero. This is not a gauge-fixing 
but a statement that the gauge field configuration is static and, consequently, 
the magnetic fields absent. d)Ao finite everywhere and zero at infinity. If we 
also use the initial ansatz for the scalar field we take for the Lagrangian: 

L = V(e - g A*) 2 S - U{#3) + \{d a A,f + (d a J)(d a a) , (23) 



5 



where Latin indices run from 1 to 3 and Greek indices run from to 3. 
Defining 

= e - gA , (24) 
regarding e as a constant and after a little algebra we take 



J U v 



(25) 



where the prime denotes the derivative with respect to the radius r and 
the last term can be absent if a{r) is a step-function. The Euler-Lagrange 
equations for the above fields <r, 6, are 



at/ 



o 



0" + ^ _ lg*(fftff)0 = . 

r 2 

The total energy of the field configuration can be written 



E = in / r dr 



1 

2?' 



(26) 
(27) 

(28) 



Now we will find the total charge. A component of the current is defined 



as: 



JaQ 



dL 


dL 




' |r« 






d(8o0) 


0(800") 




-K. 







(29) 



which in the special case of the local symmetry can be written 





1 * 

2 a 



0* 



After some algebra we take 



Qi = 2vr y r 2 drQi^B . 



Another useful relation can be obtained from equation E3 

(r 2 ey = \g\Jd)er 2 . 



(30) 
(31) 

(32) 
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By the definition of 6 and by substituting eq. ED in EH we can take an 
asymptotic relation connecting 0^, or equivalently A i0 and Qi, when r — > oo. 
The relation is: 



e i = £i -g 2 ^. (33) 

Also for large r, a is small, so as to have a localized solution of the equations 
of motion, so U((T a) = m 2 a*a and equation (J2fij) takes the form: 

a" + -a' + a( -e 2 - m 2 ) = , (34) 
r \4 / 

the solution of which has the asymptotic form: <7j oc exp (^—r^jm 2 — jS 2 ^j /r, 

with <Tj the i component of the a doublet. In order to avoid oscillatory 
solutions we demand that \e 2 < m 2 . Another useful relation can be taken 
from the asymptotic forms of a and 6 and from relations ()28j) and (|31j) : 

E = ~e • Q + 4vr y r 2 rfr[[/(o ;t a) + ^'a'] , 

where Q = (Qi,Q2,Qs) as defined above. So, if this energy is less than 
the energy of the free particles the soliton is stable against decaying to free 
particles. 

We will now try to find a solution corresponding to a quite large soliton 
in the thin-wall approximation. We will also take some general properties 
of the soliton. The more convenient method for solving the Euler-Lagrange 
equations is to regard the "matter" field a as a constant within a certain 
volume. In order to minimize the surface contribution to the total energy we 
regard the field configuration as spherical. The final scope of our work is to 
find an expression of the energy with respect to the charge, as we do in the 
global SU(2) case. We will minimize the energy firstly with respect to the 
radius of the soliton and then with respect to the field value, as is the usual 
practice in the treatment of q-balls. We will use the following approximations: 
a)The thin-wall approximation for the matter field a. b)The assumption of 
small coupling constant, g. We will not use the approximate asymptotic 
relation for the energy given above as in where the case of gauged U(l) 
q-balls was treated, but the exact one, (equation EHl)- Let the matter have a 
spatially constant value <To within a sphere of radius R. We will use eq. (31) 
so as to find 6i(r). It is easy to find that: 



(^-^f)i?sinh(^r) 



9i(r) = + . vg , , r<R (35) 
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and: 

e i {r)=e i - 9 ^, r>R, (36) 

1 /2 

where E=[(1/2)ctV] . It is easy to see that the gauge field at large dis- 
tances of the soliton origin has the form of a field generated from a spherically 
symmetric localized distribution carrying SU(2) charge, without magnetic 
fields. Eq. ED gives 

e - ^ ( 1 ^ (37 ) 

4tt \gZR - tanh(gZR) J ' v ; 

Inserting from the above equation in the relations 13*6113*71 and substituting 
the result in the exact expression of the energy we find: 

E = 7TF + W" 1 + exp(2^S)]/Q 2 - 
32[— 1 + exp(2^ J RS)]7r 2 J R 4 f/ + 32[1 + exp{2gRE)}g7i 2 R 5 m} 

{24tlR[1 + + exp{2gRE){-l + ^iffi)]} -1 . (38) 

In the above expression the first term is the energy contribution of the gauge 
field outside the q-ball. The other terms give the energy of the interior of 
the soliton. We differentiate the above expression with respect to the radius 
of the configuration and we set the result equal to zero. In order to solve the 
transcendental equation that we take, we expand it in a power series of the 
coupling constant g. The equation that we have to solve is: 

9Q2 +4^-SS = 0. (39) 



The solution to the equation should also be expanded in a power series of g: 

R 1/3\ 1/3 (2QW) 1/6 | Q 2 SW 

' ~2W ZVU Sv / f7607r(2Q 2 S 4 f/ 2 ) 1/2 ' 

Inserting the above expression into the expression of energy and expanding 
again the energy with respect to g we find 



1/2 



^ 2 (Q 2 ) 1/3 (^S) 
i + v ^ ' 



\ 2/3 



(41) 
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The limit of the above expression when (7 — > is what we expect from 
the global case. The final step is to minimize the energy with respect to the 
field S. Let S be the value of the field minimizing the energy in the global 
case. Then the value that minimizes the energy in the local case is slightly 
different from the first one, as one can expect when the coupling constant is 
small and is given by: 




d 2 fU_ 
<9£ 2 VS 2 



-r 



(42) 



where: 



e o = 



2t/(Ec 

Y 2 
^0 



(43) 



the value of the "frequency" that minimizes the energy in the g — > case. 
Let 

k = 





' d 2 








_<9E 2 




S=S 



If we substitute eq. 021 into EI] and keep only the g 2 terms of a Taylor 
expansion we take 




1 + 



3E^ /3 k(Q 2 ) 1/3 + 4 



— — 



1/3 



10 * 12 1 / VW/s 



(44) 



In conclusion, we minimized the energy with respect to the radius of the field 
configuration (eq. l4~Tj) and the field a (eq. EJ) . The energy seems to be larger 
than in the global case due to the terms depending on g 2 . There is actually 
a slight decrease in the energy with the coupling constant. This result will 
be verified numerically. But the energy also depends on the charge. This 
total charge, a/ Q 2 , may be different in the local case from the global one. 
Intuitively, we expect that this charge (particle number) is smaller in the case 
of the local symmetry due to the electrostatic-type repulsion. This result will 
also be verified by numerical analysis. 



4 Numerical results 



The Euler-Lagrange equations of the matter and the gauge field are eqs. EB1 
l2~71 The energy and the soliton charge are given from the equations l2*Hl and l3*T1 
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0.8 1 1.2 1.4 1.6 1.8 2 

Frequency 

Figure 1: The central value of the field as a function of \/~£* for three different 
values of the coupling constant with B = 0.25. 

respectively. In the global case we put g = and = e. The soliton charge 
is half the energy of the free particles with the same charge if the mass is 
unity as we have seen. 

FiguresHJ |2H2] and 0] show the behavior of the soliton as a function of \fe* . 
The square of the coupling constant takes the values 0, 0.004 and 0.04. The 
first is the case of the global symmetry. We see that when approaches the 
critical (minimum) value, y/s^, which holds for the thin-wall approximation, 
the "peak" of the soliton approaches the 2/35 value according to eq. [T3 
and the energy and the charge of the soliton increase. These all agree with 
the idea that when the frequency approaches the critical, minimum, value, 
the contribution of the spatial field derivatives to the total soliton energy get 
smaller and we can talk for a large soliton with large energy contributions 
from the time dependence and the potential. When \/s^ — > 2m the size, 
energy and charge of the soliton decrease, tending to a flat field configuration 
covering the whole space with a very small value which corresponds to free 
particles. For the above reason the ratio soliton energy per free particles 
energy approaches unity. When the coupling constant differs from zero (i.e. 
when referring to the case of local symmetry) the soliton energy and charge 
are less than in the global case. 

Figures OS] give the central value, the energy of the soliton, the energy 
of free particles with the same charge and the ratio of the two energies. The 
potential now is less "deep", as the new value of the parameter B is 0.75. 
This has some interesting consequences. The minimum value of the frequency 
(2-^/1 — 2/95) is now larger than in the case of B = 0.25. So the ratio of the 
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B=0 .25 
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Figure 2: The energy of the soliton as a function of Ve* for three different 
values of the coupling constant with B = 0.25. 




B=0 .25 
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Frequency 

Figure 3: The energy of the free particles with charge equal to the soliton 
charge as a function of Ve^ for three different values of the coupling constant 
with B = 0.25. The charge of the soliton is half the value of the energy of 
the free particles. 
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Figure 4: The soliton energy per the energy of free particles with equal charge 
as a function of \fe^ 
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Figure 5: The central value of the field as a function of Ve^ for three different 
values of the coupling constant with B = 0.75. 
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B=0 . 75 



1.7 1.75 1.8 1.85 1.9 1.95 
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Figure 6: The energy of the soliton as a function of \fe^ for three different 
values of the coupling constant with B = 0.75. 




B=0 . 75 



1.7 1.75 1.8 1.85 1.9 1.95 

Frequency 

Figure 7: The energy of the free particles with charge equal to the soliton 
charge as a function of \fe^ for three different values of the coupling constant 
with B = 0.75. The charge of the soliton is half the value of the energy of 
the free particles. 
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Figure 8: The soliton energy per the energy of free particles with equal charge 



soliton energy per the energy of the free particles is in general larger than in 
the case of B = 0.25, because the minimum frequency is the lower limit of 
the above ratio. 

In figures IUHT21 we study the behavior of the soliton when the coupling 
constant varies. When the coupling constant increases the soliton value de- 
creases. This is an expected behavior due to the repulsion between the 
different parts of the soliton and the main reason for the decrease of the par- 
ticle number (i.e. the soliton charge.) We will now explain the behavior of 
the soliton energy as a function of the coupling constant. Firstly the energy 
decreases with respect to the coupling constant due to the decrease of the 
value of the scalar field. Then, the energy seems to increase slightly due to 
the increase of the "electrostatic" energy contribution. 

5 Conclusions 

We found the energy and charge for a specific field configuration correspond- 
ing to the q-ball ansatz (eq. EJ) in a Lagrangian of a spinor doublet of two 
scalar fields with a global SU(2) symmetry. In the thin-wall approximation 
we found that, for potentials raising near the origin as |^| , then slower than 
\if\ 2 and for large \<p\ faster than \ip\ 2 , q-ball-type solitons can be observed 

at or near the minimum of the J U/\(p\ 2 quantity. In this case the soliton 
energy takes a simple form and is smaller than the energy of the free par- 
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Figure 9: The central value of the field as a function of the coupling constant 
for three different values of the frequency with B = 0.25. 
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Figure 10: The energy of the soliton as a function of coupling constant for 
three different values of the frequency with B = 0.25. 
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Figure 11: The energy of the free particles with charge equal to the soliton 
charge as a function of the coupling constant for three different values of the 
frequency with B = 0.25. The charge of the soliton is half the value of the 
energy of the free particles. 
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Figure 12: The soliton energy per the energy of free particles with equal 
charge as a function of the coupling constant with B = 0.25. 
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tides with the same charge, provided that the above quantity is less than 
the mass of the free particles. In the thick-wall approximation we used the 
full Euler-Lagrange equation for the field and after solving it numerically we 
found that for a large region in the parameter space the soliton is prevented 
from fission into free particles thanks to the energy and charge conservation. 

In the case of local SU (2) symmetry we found that a localized field config- 
uration is possible as a solution of the Euler-Lagrange equations. We found 
the energy of such a distribution as a function of charge. We minimized the 
energy with respect to the radius. Minimizing the energy with respect to the 
field value we found that this value is less than the corresponding to the global 
case due to the repulsion between the different parts of the charged soliton. 
For small values of the coupling constant and using numerical methods we 
found that a)there are localized and, thus, non-topological-soliton-type solu- 
tions to the equations of motion and b)their energy is less than the energy 
of the free particles. The comparison of the properties of the numerical solu- 
tions with different values of the parameters (coupling constant, parameters 
of the potential etc) helps us to fully understand the behavior of the soliton 
and the differences between the local and the global case. 
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